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Abstract 


We  have  developed  a  measurement  method  to  monitor  P-  and  S-waves  generated  from 
small  laboratory-scale  explosions  in  meter-sized  rock  media  at  a  series  of  stations  as 
well  as  invented  a  unique  device  to  drill  spherical  cavities  with  the  diameter  up  to  10 
centimeters  in  rock. 

We  have  also  conducted  laboratory  experiments  in  which  spherical /cylindrical 
charges  in  0.2  to  1.9  gram  mass  range  were  detonated  at  the  center  of  spherical/ cylindrical 
cavities  with  the  radius  from  0.6  to  1.5  cm  in  Bedford  limestone  and  the  waves  gen¬ 
erated  by  explosions  were  recorded  in  the  near  field  within  250  mm.  The  records 
and  post-explosion  studies  clearly  demonstrate  the  P-waves  generated  from  both  ex¬ 
plosions  in  cavities  and  crack  propagation  in  rock.  The  peak  particle  displacement 
induced  by  P-waves  ranges  from  9  to  5  /im  at  150  to  210  mm  from  source  for  the 
spherical  cavity  experiment  with  1.9  gram  explosives,  and  8  to  4  /zm  at  150  to  220  mm 
from  source  for  spherical  experiment  with  0.8  gram  explosives,  35  to  10  /zm  at  150 
to  240  mm  from  source  for  the  tamped  cylindrical  cavity  experiment  with  1.2  gram 
explosives  and  7  to  2  ^m  at  150  to  190  mm  from  source  for  the  decoupled  cylindrical 
cavity  experiment  with  0.24  gram  explosives.  The  step-like  pressure  model  can  be 
used  to  explain  the  experimental  results  on  the  peak  P-wave  induced  particle  dis¬ 
placement  versus  propagation  distance  from  the  explosions  in  spherical  cavities  but 
not  in  cylindrical  cavities.  This  results  demonstrates  the  cavity  geometrical  effects  on 
wave  pattern.  The  experimental  results  from  both  spherical  and  cylindrical  cavities 
show  that  P-wave  profiles  are  similar  with  each  other.  The  experimental  results  verify 
that  S-waves  are  generated  from  both  shock  wave-to-S  conversion  at  cavity  walls  and 
crack  propagation  in  rocks.  The  amplitude  of  S-waves  from  shock  wave-to-S  conver¬ 
sion  is  from  10  to  27  /zm  at  150  to  180  mm  from  source  for  the  spherical  experiment 
with  1.9  gram  explosives,  8  to  21  /zm  at  150  to  180  mm  from  source  for  the  spherical 
experiment  with  0.8  gram  explosives.  The  amplitude  of  S-waves  generated  by  crack 
propagation  is  from  9  to  19  /zm  at  150  to  180  mm  from  source  for  the  spherical  ex- 
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periment  with  1.9  gram  explosives,  8  to  19  /xm  at  150  to  180  mm  from  source  for  the 
spherical  experiment  with  0.8  gram  explosives.  Based  on  the  experimental  data,  the 
ratio  of  S-wave  to  P-wave  amplitude  is  approximately  about  2. 

We  have  analyzed  P  and  S-waves  generated  by  explosions  and  crack  propagation 
in  rock.  A  theoretical  model  is  proposed  to  explain  the  differences  of  P-wave  induced 
displacement  history  between  the  observed  waveforms  and  that  predicted  by  a  step 
pressure  source.  The  predications  of  this  model  fit  the  observation  qualitatively.  The 
present  results  demonstrate  the  important  role  of  rock  cracking  in  P  and  S-wave 
generations  from  overdriven  cavities. 

The  results  from  this  work  have  demonstrated  that  the  seismic  source  for  under¬ 
ground  explosions  in  cavities  is  much  more  complex  than  seismic  theoretical  models 
used  to  describe  underground  explosions.  The  rich  phenomena  associated  with  waves 
generated  from  explosive  product  driven  crack  extending  from  cavities  suggest  that 
this  research  direction  may  lead  to  development  of  the  criteria  that  can  be  used 
to  discriminate  underground  explosions  from  earthquakes  and  possible  discriminate 
chemical  from  nuclear  explosion  sources. 


Chapter  1  Introduction 


In  the  past  30  years,  extensive  research  has  been  conducted  on  the  discrimination 
of  underground  explosions  from  earthquakes.  As  a  result,  Murphy  (1996)  recently 
concluded  that  the  seismology  of  earthquakes  and  underground  explosions  is  at  a  state 
that  it  is  possible  to  identify  tamped  explosions  with  a  yield  over  1  kiloton  (  mj  > 
4)  on  a  world-wide  basis.  These  explosions  can  be  distinguished  from  earthquakes 
by  the  global  seismic  stations  using  the  initial  P-wave  motion  and  the  M,  to  mb 
ratio.  In  fact,  the  discrimination  of  the  tamped  explosions  from  earthquakes  depends 
mainly  on  the  collection  of  seismic  data  from  previous  underground  nuclear  explosions. 
However,  many  basic  problems  related  to  tamped  explosions  still  remain  unsolved  ( 
Masse,  1981;  Murphy,  1996  ),  e.g.,  the  generation  of  long-period  Rayleigh  waves  (>  5 
seconds),  SH- waves  and  Love- waves  etc.  But  now,  the  most  difficult  problems  are 
to  distinguish  earthquakes  from  decoupled  explosions  with  complex  geometries  and 
chemical  explosions  from  small  nuclear  explosions. 

The  characteristics  for  decoupled  explosions  (Murphy,  1996)  are:  (1)  my  is  too  low 
to  use  the  global  monitoring  system  ( in  general,  mj,  <  3);  (2)  the  amplitudes  of  surface 
waves  generated  from  such  explosions  are  much  lower  than  from  tamped  explosions; 
(3)  the  excitation  of  direct  SV-waves  and  possibly  SH-waves  makes  the  problem  more 
difficult.  Although  there  are  many  published  papers  related  to  decoupled  explosions 
(Glenn  et  al.,  1985, 1993, 1994, 1996;  Murphy  et  al.  1996;  Murphy,  1996;  Sykes,  1996; 
Langston,  1983  etc.),  it  still  is  unclear  what  really  controls  the  generation  of  S-waves 
from  decoupled  explosions  (Murphy,  1996).  We  conclude  that  the  mechanisms  of 
S-wave  generation  from  explosions,  especially  from  decoupled  explosions  in  cavity  is 
potentially  an  important  discrimination  criterion. 

The  suggested  possible  sources  for  S-waves  from  explosions  include: 

1.  dynamic  fracture  of  rocks  near  explosions; 

2.  the  asymmetric  plastic  deformation  near  the  explosion  due  to  the  asymmetric 
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geometry  and  anisotropy  of  the  rocks; 

3.  tectonic  release  of  shear  energy  triggered  by  shock  waves  from  the  explosions; 

4.  the  conversion  of  P-waves  to  S-waves  at  the  free  surface  of  the  Earth; 

These  mechanisms  must  be  studied  separately  using  different  methods. 

The  objective  of  this  work  is  to  investigate  P-  and  S-wave  generations  from  ex¬ 
plosions  in  cavities,  dynamic  fracture  in  rock  and  the  conversion  of  shock  induced  P- 
and  S-waves  at  the  surfaces  of  the  cavities  experimentally.  If  the  shock  wave  pres¬ 
sure  in  cavity  is  not  too  high,  the  rocks  near  explosions  will  fail  via  brittle  fracture 
mechanisms.  Radial  fractures  may  be  initiated  from  the  inner  surface  of  the  cavity 
because  the  tangential(hoop)  stresses  always  achieve  the  maximum  value  at  the  sur¬ 
face.  Radial  fractures  generally  can  be  further  driven  by  the  high-pressure  explosive 
products  (Coursen,  1985).  The  propagation  of  radial  cracks  may  result  in  P-  and 
S-wave  generations.  The  reflection  of  the  shock  waves  at  cavity  surfaces  will  also 
induce  S-waves. 

In  order  to  investigate  the  details  of  the  P-  and  S-wave  generations  from  explosions 
in  cavities  in  rock,  a  measurement  method  that  has  the  capability  to  measure  P-  and 
S-waves  generated  from  explosions  in  meter-size  samples  and  a  device  that  can  form 
spherical  cavities  in  rock  have  to  be  developed. 

For  small-scale  laboratory  experiments,  the  wave  proflles  have  much  higher  fre¬ 
quencies  than  the  usual  seismic  waves.  Therefore,  conventional  seismic  recording 
systems  can  not  be  used  to  measure  those  wave  proflles.  One  candidate  method  de¬ 
veloped  for  super-high  strain  rate  experiments  we  considered  (Kim  and  Clifton,  1977) 
is  difficult  to  employ  for  the  present  experiments  because  it  requires  special  treatment 
of  sample  surfaces.  Magnetic  velocity  gauges  have  long  been  used  to  study  shock  wave 
decays  in  rocks,  but  this  method  requires  embedding  the  gauges  in  the  rock  and  thus 
contains  considerable  machining  of  rocks.  This  also  adds  some  artificial  internal  in¬ 
terfaces  which  will  generate  S-waves  and  make  the  S-wave  train  more  complicated. 
Based  on  the  analysis  of  the  interaction  between  P-  and  S-waves  and  free-surfaces, 
we  have  developed  a  method  to  measure  P-  and  S-wave  profiles  on  the  free  surfaces 
of  rocks.  The  main  idea  of  the  method  is  to  use  the  different  characteristics  of  the 
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interactions  between  P-  and  S-waves  and  the  free  surface.  In  order  to  make  a  spherical 
cavity  in  rock,  a  special  drilling  tool  has  been  developed  and  used  in  the  experiments 
successfully. 

In  this  report,  the  details  of  the  method  including  measuremental  method  and 
spherical  cavity  drilling  method  are  described.  Experimental  data  on  spherical  and 
cylindrical  cavities  and  the  analysis  on  P-  and  S-wave  profile  and  amplitudes  are 
presented. 
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Chapter  2  Measurement  Method 


2.1  Introduction 

When  P-  and  S- waves  arrive  at  a  free  surface,  the  interactions  between  the  P-  and  S- 
waves  and  the  free  surface  are  very  different.  The  characteristics  of  these  interactions 
provide  the  opportunity  to  measure  the  incident  P-  and  S-wave  amplitudes  from  the 
strain  measurements  on  the  free  surfaces  of  rocks.  The  method  developed  is  to  use 
two  strain  gauges  to  measure  the  strains  along  two  perpendicular  directions  at  one 
point  on  the  free  surface  of  rock.  The  incident  amplitude  of  P-  and  S-waves  can  be 
obtained  from  these  strain  measurements. 

The  basic  assumptions  made  for  the  method  are: 

1.  The  waves  in  the  rocks  can  be  approximately  treated  locally  as  plane  waves 
because  the  dimensions  of  strain  gauges  used  in  experiments  are  only  3x3  mm 
and  gauges  are  located  at  least  150  mm  away  from  the  explosive  sources. 

2.  The  deformation  of  the  rocks  near  the  free  surface  is  elastic. 

The  following  is  the  details  of  the  derivation  of  the  relations  between  the  signals 
from  strain  gauges  and  P-  and  S-waves  in  rock. 

2.2  Sensors 

A  group  of  strain  gauges  are  employed  to  monitor  the  wave  profiles  at  different  loca¬ 
tions  along  two  directions  on  the  free-surface  of  rock  samples.  The  directions  of  the 
gauges  are  shown  in  figure  2.1.  The  gauges  and  amplifiers  are  powered  by  batteries 
in  order  to  achieve  a  high  signal  to  noise  ratio.  A  mercury  relay  is  used  to  calibrate 
the  system.  The  relation  between  strains,  e,  and  the  voltage  change,  AV,  is 


where  V©  is  the  initial  voltage  supplied  on  the  gauges.  /  and  g  are  the  gauge  factor 
and  gain  of  the  amplifiers,  respectively. 


2.3  Data  Reduction  Method 

Strain  gauges  give  the  elongation  along  the  gauge  directions.  When  they  are  attached 
on  a  free  surface,  they  measure  the  elongation  of  the  free-surface  after  wave  reflections. 
The  strains  recorded  by  the  gauges  include  the  contributions  from  incident  P-  and 
S-waves  and  the  reflected  P-  and  S-waves.  In  order  to  get  the  relation  between  the 
strains  given  by  the  gauges  and  the  incident  P-  and  S-wave  amplitudes,  we  first 
applied  the  well-known  formulas  describing  the  interactions  between  the  incident  P- 
and  S-waves  and  the  free  surface. 


2.3.1  P-wave  reflection  at  free  surfaces 


The  displacement  reflection  coefficients  for  incident  P-waves  at  free  surfaces  (Aki  and 
Richards,  1980)  are 


PP  = 


B-A 
B  +  A' 


PS  = 


2^  sin(20)  cos(2j) 
A  +  B 


(2.2) 

(2.3) 


where  PP  and  PS  are  the  reflection  coefficients  for  P-  and  SV-wave  displacements 
due  to  the  incident  P-waves,  respectively,  a  and  P  are  P-  and  S-wave  velocities, 
respectively.  A  and  B  are 

A  =  cos^(2j),  (2.4) 

B  =  (— )^sin(2j)sin(20),  (2.5) 

a 

where  6  and  j  (see  2.1)  are  the  P-wave  incident  angle  and  S-wave  reflection  angle, 
respectively. 


.  .  ^ 

sinj  =  —  sinw, 

a 

cosj  =  (1  -  (-)^sin^  j)2. 
a 


(2.6) 

(2.7) 
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The  resultant  displacements  of  the  particles  on  the  free  surfaces  after  the  reflection 


are 


+  PP)  sin0  +  PS  cos;]  =  ifparUp, 

(2.8) 

=  Up[(l  -  PP)  cose  +  PS  sin;]  =  HperU^, 

(2.9) 

where  is  the  particle  displacement  of  the  incident  P-wave,  and  indicate 

the  resultant  particle  displacements  along  direction  2  and  the  direction  that  is  per¬ 
pendicular  to  the  free  surface  after  reflection  as  shown  in  figure  2.1,  respectively.  The 

PpaT  and  Hpof  are 

Hpar  =  (1  +  PP)  sin^  -t-  PS  cos;, 

(2.10) 

Hper  =  (1  -  PP)  cos6  +  PS  sin;. 

(2.11) 

Substituting  PP  with  Eq.(2.2)  and  PS  with  Eq.(2.3),  Eqs.(2.10)  and  (2.11)  can 

be  rewritten  as 

^  2cos0sin(2;) 

A^B  ' 

(2.12) 

2cos0cos(2;) 

A  +  B  • 

(2.13) 

2.3.2  SV-wave  reflection  at  free  siurfaces 

For  incident  SV-waves,  the  displacement  reflection  coefficients  for  P- 

and  SV-waves 

(Aki  and  Richards,  1980)  are 

A, B,' 

(2.14) 

ss=^‘~^\ 

A,  -f  Bg 

(2.15) 

where  SS  and  SP  are  the  reflection  coefficients  of  SV-  and  P-wave  displacements  due 

to  the  incident  SV-waves. 

Ag  =  cos^(2^). 

(2.16) 

B,  =  (-)^sin(2^)sin(2;), 
a 

(2.17) 
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where  6  and  j  are  the  SV-wave  incident  angle  and  P-wave  reflected  angle,  respectively, 
and 


cosj  =  (1  - 

(2.18) 

The  resultant  displacements  of  the  particles  on  the  free  surfaces 

after  the  reflection 

are 

»^par  =  +  *55)  COS^  +  SPsinj]  -  Gparui„, 

(2.19) 

“per  =  “L[(‘S'5  -  1)  sine  -  5P cos i]  =  Gperui^, 

(2.20) 

where  u*",.  and  are  the  resultant  particle  displacements  on  the  free  surface  after 

reflection  along  direction  2  as  shown  in  figure  2.1  and  the  direction  that  is  perpendic¬ 
ular  to  the  free  surface,  respectively.  is  the  particle  displacement  of  the  incident 

S\^— wave.  The  Gtp^^  and  ^^por  are 

Gpar  =  (1  +  SS)  COS  6 SP  sin  j, 

(2.21) 

Gper  —  {SS  -  1)  sin  9  -  SP  cosj. 

(2.22) 

Substituting  SP  with  Eq.(2.19)  and  SS  with  Eq.(2.15),  Eqs.(2.21)  and  (2.22)  are 

rewritten  as 

2  cos(20)  cos  6 
“  As +  B,  ’ 

(2.23) 

^  cosisin(20) 

As  +  Bs  ‘ 

(2.24) 

2.3.3  Strains  induced  by  incident  P-waves 

1.  Strains  measured  along  1  direction 

Because  is  perpendicular  to  the  free  surface  and  along  the  symmetrical  axis 
direction  of  the  waves,  the  strain  due  to  along  the  1  direction  is 

eT  =  (2-25) 

ro 

where  ef®’'  is  the  strain  along  direction  1  induced  by  is  the  distance 

from  the  center  of  the  cavity  to  the  free  surface  at  0  =  0. 
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Because  does  not  result  in  any  strains  in  the  gauges  along  direction  1  at  any 
time,  the  total  strain  induced  by  the  incident  P-waves  is 

ef  =  (2.26) 

To 

where  Hi  =  Hper- 

2.  Strains  measured  along  2  direction 

Because  and  all  have  contributions  to  the  strain  measured  by  gauges 
along  the  2  direction,  we  need  to  consider  the  resultant  displacements. 

The  length  of  the  gauge  after  the  reflection,  V,,  is 

r.  =  (r;  +  (2.27) 

where 

59=—,  (2.28) 

Tn 

and  It  is  the  initial  length  of  strain  gauges.  Vn  is  the  distance  from  the  center  of 
the  cavity  to  the  position  of  the  gauge  after  reflections 


r„  =  r  +  ucos(77  —  9), 


(2.29) 


where  r  is  the  distance  from  the  center  of  the  cavity  to  the  gauge  before  P-wave 
reflection,  u  is  the  resultant  displacement  of  the  point  at  9  on  the  free  surface 
and  T]  is  the  angle  between  u  and  u^,  they  are 


u  =  u: 


j2cos(0) 
^  A  +  B' 


(2.30) 


V  =  2i. 


(2.31) 


To  a  first-order  approximation,  is 

((^)"  +  +  (2.32) 
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where 


2  cos(0)  cos(77  —  9) 

~  A  +  B 

(2.33) 

If  X  ^  1,  we  have 

(l  +  x)*  ~  1  +  ^. 

(2.34) 

Because  Up/r  <  1, 

after  using  the  approximation  above,  is 

i:  «  (r(l  + 

+  u[(W{e)(l  -  ^^)  +  tan 

(2.35) 

Then  the  strain  is 

^2—  7  ? 

(2.36) 

and  therefore 

e^2  =  ^2^, 

»’o 

(2.37) 

where 

„  (cos  9(W'(l-af«)  + tans®) 

2  (1  +  i^) 

(2.38) 

and  ro  =  r  cos  6. 


2.3.4  Strains  induced  by  incident  S- waves 

Using  the  same  method,  all  the  strains  induced  by  incident  SV-waves  can  be  obtained. 
The  final  expressions  are  listed  as  following: 

1.  Strains  measured  by  gauges  along  1  direction  are 

er  =  C?i— ,  (2.39) 

To 

where 

~  ^per-  (2.40) 

2.  The  strains  given  by  the  gauges  along  direction  2  can  be  found  using  the  same 
method  as  for  the  incident  P-wave.  The  resultant  strains  induced  by  the  incident 
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SV-waves  along  direction  2  are 


-$V  _  '"»v 

^2  ~  ^2 - j 

^■o 


where 


G2  = 


(cos^(W,(l 


0  +  tan^^) 


(l^ta^) 


W,= 


2cos(77-  ^)cos^.  2/onx  .  ..„:„2mi 


Ai  +  Bg 

tanrj  = 


(cos^(2^)  +  4(— cos^  j  sin^  0)  * , 
a 

a  cos  6 


P cosj tan(2^) 


(2.41) 


(2.42) 

(2.43) 

(2.44) 


2.3.5  Amplitudes  of  incident  P-  and  S-waves 

From  the  expressions  given  above,  the  displacements  of  incident  P-  and  SV-waves  can 

be  obtained  through  the  strains  given  by  the  gauges  along  the  two  directions. 

1.  The  displacement  of  the  incident  P-waves  is 

=  ’■<’§;  = 

2.  The  displacement  of  the  incident  SV-waves  is 

<  =  ro^  =  (2.46) 

2.4  Characteristics  of  the  Strains  along  the  Two  Directions 

2.4.1  Strains  induced  by  incident  P-waves 

1.  The  non-dimensional  constant,  Hi  in  Eq.(2.45)  is  not  sensitive  to  the  variation  in 
9,  e.g.,  the  value  of  Hi  is  about  2.0  for  Bedford  limestone.  The  non-dimensional 
constant,  H2,  in  Eq.(2.45)  changes  rapidly  with  $,  it  varies  from  about  2  to  about 
-0.8  when  incident  angle  changes  from  0  to  55°  for  Bedford  limestone  as  shown 
in  figure  2.2. 


10 


61  cm 


Strain  gauges  (7,  8) 


(A)  Layout  of  rock  sample,  strain  gauges  and  explosive  charge 


(B)  Sketch  of  polarization  directions  of  strain  gauges 
Figure  2.1:  Sample  and  strain  gauge  layout  for  spherical  cavity  experiments 
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2.  Prom  figure  2.2,  the  strains  induced  by  P- waves  along  direction  1  are  always 
positive,  but  the  strains  along  direction  2  are  positive  when  6  is  less  than  47° 
and  negative  when  B  is  larger  than  47°  (  for  Bedford  limestone).  This  polarity 
change  is  controlled  by  the  ratio  of  the  projection  of  P-wave  displacement  cilong 
direction  1  to  that  along  the  direction  that  is  perpendicular  to  the  free  surface. 

3.  The  characteristics  of  the  strain  induced  by  incident  P-waves  have  been  verified 
by  the  experimental  results  as  shown  in  figures  3.5,  3.6,  4.7  and  4.8.  Prom  the 
experimental  results,  it  is  seen  that  all  the  strains  along  direction  1  (labeled 
with  even  gauge  number)  are  positive  at  all  the  stations  and  the  strains  along 
direction  2  (labeled  with  odd  gauge  number)  are  positive  when  the  incident  angle 
is  smaller  than  the  critical  angle  and  negative  when  the  incident  angle  is  larger 
than  the  critical  angle. 

2.4.2  Strains  induced  by  incident  S-waves 

1.  The  relation  between  G\  {G-i)  and  the  SV-wave  incident  angle  for  Bedford  lime¬ 
stone  is  given  in  figure  2.2.  From  the  calculated  results,  the  gauges  along  direc¬ 
tion  1  are  not  sensitive  to  incident  SV-wave;  however,  the  gauges  along  direction 
2  are  very  sensitive  to  incident  SV-wave. 

2.  The  measurement  of  SV-waves  are  limited  by  the  Rayleigh  surface  wave  gener¬ 
ation.  For  Bedford  limestone,  the  Rayleigh  wave  is  generated  when  6  is  larger 
than  about  35°  (  determined  from  Snell’s  law). 

3.  The  polarities  of  the  strains  along  direction  2  are  negative,  and  the  polarities 
of  the  strains  along  direction  1  are  determined  by  the  direction  of  the  particle 
motion  (  the  calculation  is  made  assuming  that  the  motion  direction  is  toward 
the  increasing  direction  of  the  SV-wave  incident  angle). 

4.  The  characteristics  of  the  strain  induced  by  incident  S-waves  have  been  verified 
by  the  experimental  results  as  shown  in  figures  3.5,  3.6,  4.7  and  4.8.  From  the 
experimental  results,  it  is  seen  that  all  the  strains  along  direction  1  (labeled  with 
even  gauge  number)  are  generally  very  small  compared  to  the  correspond  strain 
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Incident  angle  (degree) 

Figure  2.2:  H  and  G  versus  incident  angle 

amplitudes  along  direction  2  at  all  the  stations  and  the  strains  along  direction  2 
(labeled  with  odd  gauge  number)  are  always  negative  at  all  the  stations. 

2.5  Conclusions 

A  measuremental  method  that  can  be  used  to  monitor  P-  and  S- waves  generated  from 
explosions  has  been  developed.  The  relations  between  incident  P  and  S-wave  particle 
displacement  and  the  strains  recorded  along  the  two  directions  have  been  derived  as 
given  in  Eqs.(2.45)  and  (2.46)  and  the  experimental  set-up  is  shown  in  figures  2.1  and 
3.3.  Experimental  results  have  verified  the  P-  and  S-wave  characteristics  predicted 
from  the  derived  relations  as  shown  in  figure  2.2. 
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Chapter  3  Wave  Generation  Prom  Explosions  in 
Spherical  Cavities  in  Rock 


3.1  Introduction 

Nearly  all  the  seismic  models  for  underground  explosions  assumed  that  the  source 
is  spherical  (Murphy,  1996).  Therefore,  experimental  study  on  the  P  and  S-waves 
generated  from  the  explosions  in  spherical  cavities  has  direct  applications  to  the 
understand  of  the  waves  from  underground  explosions.  The  objectives  of  this  part  is 
to  study  the  P  and  S-waves  generated  from  explosions  in  spherical  cavities  and  also 
from  crack  propagation  driven  by  shock  waves  and  the  explosive  products. 

3.2  Experiments  and  Results 

3.2.1  Spherical  cavity  drill 

In  order  to  drill  spherical  cavities  in  rock,  a  special  tool  has  been  developed.  The 
device  is  schematically  shown  in  figure  3.1.  Figure  3.2  is  the  assembly  drawing  of  the 
drill.  The  principal  of  the  drill  is  to  use  one  of  the  two  shanks,  Shank  A,  to  transfer 
the  force  from  the  machine  to  the  other  shank.  Shank  B,  that  rotates  the  blade  and  at 
same  time  push  the  blade  to  cut  rock.  During  drilling  process,  the  positioning  plate 
is  used  to  keep  the  drill  from  vibration  induced  by  the  asymmetrical  force  acting  on 
the  shank  from  blade.  Pressurized  air  is  used  to  cool  the  blade  and  also  blow  out 
the  rock  dust.  Principally,  the  device  can  be  used  to  drill  spherical  cavities  with 
the  diameter  up  to  10  centimeters  in  rocks  like  limestone.  The  procedures  to  form  a 
spherical  cavity  inside  rock  are 

1.  We  make  a  cylindrical  hole  with  1  cm  diameter  using  a  drill.  The  depth  of  the 
hole  will  be  the  deepest  point  of  the  spherical  cavity. 

2.  After  drilling,  a  small  size  blade  with  1.5  cm  diameter  is  put  into  the  cylindrical 
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hole,  the  position  plate  is  screwed  to  the  shank  to  fix  the  depth  of  the  spherical 
cavity  and  the  center  line  in  the  rock  target. 

3.  During  drilling,  pressurized  air  is  used  to  blow  out  rock  dust  and  also  cool  the 
blade.  The  air  flow  rate  is  roughly  about  2  1/minute.  The  rotation  speed  is 
about  240  r.p.m.  and  the  feed  rate  is  about  1  mm/minute. 

4.  Then,  a  large  blade  is  used  to  enlarge  the  cavity. 

The  device  and  procedure  have  been  used  to  make  spherical  cavities  in  limestone 
successfully. 
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Figure  3.2:  Assembly  drawing  of  spherical  cavity  drill 
Blade  material:  heat  treated  tool  steel 
Shank  material:  steel 
Dimension  imit:  in 
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3.2.2  Experimental  details 


The  rock  sample  (Bedford  limestone)  is  assembled  using  two  blocks  as  shown  in  figure 
2.1.  The  rock  sample  with  strain  gauges  is  placed  inside  a  tank  pressurized  to  10  bar. 
The  type  of  strain  gauge  used  in  this  work  is  CEA-00-062UT-120  from  Measurements 
Group.  The  dimensions  of  the  gauges  are  3x3  mm.  The  initial  resistance  of  gauges 
is  120  ft  ±0.4%  and  the  gauge  factor  is  2.090  ±  0.5%.  The  procedure  to  attach  the 
gauges  on  the  surface  is:  (1)  The  surface  is  polished  with  #240  sandpaper  ;  (2)  the 
surface  is  cleaned  using  acetone;  (3)  the  gauges  are  attached  on  the  surface  using  910 
adhesive  (Permabond  International);  (4)  after  the  adhesive  is  dried,  a  thin  layer  of 
the  epoxy  (  5-minute  epoxy,  ITWDevcon)  is  applied  on  the  gauge  surface  in  order  to 
protect  them.  The  voltage  on  the  gauge  is  typical  6  V  (  before  each  experiment,  the 
voltage  on  each  gauge  is  measured).  The  amplifiers  used  in  this  work  are  specially 
designed  to  satisfy  the  requirements  of  high  gain,  wide  bandwidth  and  low  noise. 

The  explosives  used  in  this  work  are  Deta  sheet  PETN  (DuPont)  and  Pellet  pressed 
PETN  (Reynolds  Industries  System  Incorporated),  the  initial  density  is  1.5  gram/cm^ 
for  Deta  sheet  PETN  and  1.4  gram/cm^  for  Pellet  pressed  PETN.  The  size  of  the 
detonator(Reynolds  Industries  System  Incorporated)  is  (p2xS  mm. 

We  have  fired  several  decoupled  explosion  tests  in  spherical  cavities.  Due  to  inher¬ 
ent  difficulty  to  properly  detonate  small  charge,  we  put  much  effort  in  completing  a 
very  few,  highly  successful  experiments.  The  experimental  setup  is  sketched  in  figures 
2.1  and  3.3.  Figure  3.4  schematically  shows  the  spherical  explosive  source.  Table  3.1 
lists  the  parameters  for  6  shots  and  the  equilibrium  pressure,  p,  in  cavities  estimated 
using  the  expression(Murphey,  1961) 


P  = 


(7  - 1)-£^ 

y  ’ 


(3.1) 


where  7  is  the  ratio  of  specific  heats  of  explosive  products  in  air,  E  is  the  total  energy 
of  explosives  used,  V  is  cavity  volume.  In  the  calculation,  7  is  1.4  (Baker,  1973)  and 
the  energy  released  from  PETN  is  1255  J/g  (Dobratz,  1974). 
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Table  3.1:  Data  for  spherical  cavity  experiments 


Shot 

number 

Explosive 
mass  (gram.) 

P  (bar) 

SW4 

1.9 

675 

SW5* 

0.4 

142 

SW6* 

0.6 

213 

SW7’ 

0.7 

248 

SW8  * 

0.8 

284 

SW9 

0.8 

284 

*  -  indicates  that  the  explosives  did  not  properly  detonate. 


Strain  gage  cabling 

Out-going  tube  for  gage  cables 
Relief  valve 
Compressed  gas 

Rock  support  ( top ) 

Gas-in  tube 
Tank  body 

Water  ( under  pressure) 


Joint 

Rock  sample 
Plastic  film 
Pressure  gage 

Detonator  leads 

Rock  support  ( base) 

Out-going  tube  for  detonater  leads 
Water-intake  tube 

Water-outtake  tube 
Seal  nut 


Figure  3.3:  Mechanical  configuration  of  the  test  chamber. 
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Figure  3.4;  Spherical  explosive  source  arrangement. 

3.2.3  Experimental  records 

Using  the  experimental  system  given  above  and  the  measuremental  method  described 
in  Chapter  2,  the  waves  generated  from  the  explosions  in  spherical  cavity  were 
recorded.  Experimental  records  along  the  two  directions  are  shown  in  figure  3.5 
for  SW4  and  figure  3.6  for  SW9.  The  comparisons  between  the  records  at  different 
locations  are  shown  in  figures  3.7  and  3.8. 

3.2.4  Crack  pattern  &om  the  recovered  samples 

(a)  SW4 

The  recovered  experimental  sample  shows  that  there  are  four  large  cracks(~  10 
cm)  and  many  small  cracks  (~  5  mm)  induced  by  the  explosion(Fig.  3.11).  The 
length  of  the  cracks  with  the  explosion  product  trace  ranges  from  50  mm  to  80  mm. 

(b)  Shot  SW9 

From  the  recovered  sample,  cracks  were  not  seen  from  the  free  surfaces.  After  the 
sample  was  cut  along  the  axis,  many  small  cracks  near  the  cavity  existed  and  also 
two  very  long  cracks  are  formed.  The  longest  cracks  are  about  50  mm  in  length. 
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Figure  3.6:  Strain  signals  from 


Figure  3.7:  Strain  signals  at  the  same  location  from  the  gauges  in  SW4  and  SW9. 

Tsl  and  ATsl:  the  arrival  and  duration  of  the  shear  waves  induced  by  microcracking  (Eqs.(3.83)  and 
(3.84)) 

Ts2  and  ATs2:  the  arrival  and  duration  of  the  shear  waves  induced  by  shock  wave-to-S  conversion  at  the 
cavity  wall  (Eq.(3.82)) 

Ts3e  and  Ts3j:  the  possible  earliest  and  last  arrival  of  the  shear  waves  induced  by  macrocracking  (Eqs.(3.86) 
and  (3.87)) 


1\ 


Displacement  (  pm)  Strain  (P0 


3.6cm 


30.5cm 


Figure  3.9:  Schematic  of  the  large  cracks  in  relation  to  the  spherical  cavity. 

3.3  Analysis  of  the  P- waves  Induced  by  Explosions  and  Crack¬ 

ing 

3.3.1  P-wave  velocity 

From  the  experimental  records,  the  P-wave  velocity  is  determined  to  be  4.22  km/s 
under  the  experimental  conditions  as  shown  in  figure  3.10  (the  data  from  SW9  was 
not  used  because  of  the  noise  at  the  beginning  of  the  signals).  The  P-wave  velocity 
determined  here  is  somewhat  lower  than  4.9  km/s  measured  ultrasonically  by  Rubin 
and  Ahrens  (1991)  on  a  different  sample.  However,  this  velocity  is  still  much  higher 
than  the  bulk  velocity  (3.68  km/s).  This  means  that  the  deformation  induced  by  the 
P-wave  is  still  in  the  elastic  deformation  regime  over  most  of  its  propagation  path. 
The  lower  P-wave  velocity  could  also  be  the  result  of  the  nonlinear  elastic  behavior 
of  the  rock  because  elastic  wave  velocity  decreases  with  amplitude  in  rocks. 

3.3.2  Particle  displacement  induced  by  P- waves 

Using  the  relations  between  strains  and  displacements  given  in  Chapter  2,  the  P-wave 
particle  displacements  are  calculated  and  shown  in  figure  3.11  for  SW4  and  figure  3.12 


23 


Figure  3,10:  P-wave  arrival  versus  radius. 

for  SW9.  From  figures  3.11  and  3.12,  the  maximum  P-wave  particle  displacements, 
Umox.j  at  different  distances  from  the  cavity  center  can  be  obtained  and  are  shown 
in  figure  3.13  in  which  Cowboy  data  (Murphey,  1961)  are  plotted.  The  experimental 
data  is  fit  by 


Umax.  =  1-32  X  lOV-^  '*®, 

for  SW^ 

(3.2) 

Umax.  =  1.93  X  lOV-^  ®®, 

for  5IF9 

(3.3) 

where  Umax,  is  in  /xm,  r  in  mm  is  the  distance  from  the  cavity  center.  In  Mur- 
phey’s(1961)  paper,  he  found  that  the  field  results  have  a  dependence  of  the  peak 
P-wave  particle  displacement  on  the  distance  from  the  cavity  center  that  is  propor¬ 
tional  to  r“^  ®.  Thus,  the  previous  study  is  very  comparable  with  the  results  from 
this  work. 


3.3.3  Comparison  between  experimental  data  and  theoretical  predictions 
Displacement  induced  by  pulse-spherical  waves 

The  analytic  solution  for  a  pulsed  spherical  wave  generated  from  an  explosion  in  a 
spherical  cavity  (Timoshenk  and  Goodier,  1970)  is  obtained  by  solving 


da  2.  -  .dV  ^^dV, 

-  +  -K  -  a.)  =  p(—  +  V— ), 


(3.4) 
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Figure  3.13:  Scaled  maximum  P-wave  particle  displacement  versus  propagation  distance 


u  —  peak  P-wave  particle  displacement  in  mm 
w  —  mass  of  explosives  in  gram 
r  —  distance  from  the  center  of  cavity  in  mm 

For  Cowboy  experiments, 

u/w^/^  =  0.55  for  9.15  meter  cavities  in  which  the  mass  of  the  explosives  detonated 

ranged  from  200  to  1000  lbs, 

u/w^/^  =  7.7  for  3.66  meter  cavities  in  which  the  mass  of  the  explosives  detonated 

was  from  20  to  900  lbs. 

For  present  data, 

u/wl/3  ^  7  g  foj.  SW4 

u/wV3  ^  23.4  for  SW9 
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dp  1  d{pr^V) 
dt^  dr 


=  0, 


(3.5) 


where  p  is  density,  V  is  particle  velocity,  r  is  the  distance  from  center  of  the  source. 


After  using  Hooke’s  law,  the  wave  equation  is 


d^u  2du  2u  _  1  d^u 
dr^  ^  T  dr  C7|  dt^  ’ 


(3.6) 


where  u  is  particle  displacement  and  Cg  is  P-wave  velocity. 

If  pressure  at  the  cavity  boundary,  r=a,  is  p(t),  the  displacement  (Timoshenk  and 
Goodier,  1970)  can  be  calculated  using 


u  = 


C.r  r’’ 


(3.7) 


where  f  is  the  solution  of 


/"(«)  + 27/(0 +  ^/(t)  =  ^. 

a  p 


(3.8) 


where  7  is 


l  —  2vCc 
\  —  u  a  ' 


(3.9) 


where  u  is  Possion’s  ratio. 

The  particular  solution  of  Eq.(3.8)  for  an  outgoing  wave  from  a  spherical  source 
(Timoshenk  and  Goodier,  1970)  is 


m  =  rp(T)(e-"-")  -  e«"l'-">)<iT,  (3.10) 

a  —  p  p  Jo 

where 

0^0  =  7(-l  +  is),  (3.11) 

do  =  7(-l  -  is),  (3.12) 

and  s=(l  —  2v)~^'^. 

If  a  step-like  pressure  history,  i.e.,  p(t)=po  H(t),  is  assumed,  the  particular  solution 
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Figure  3.14:  Typical  displacement  history  from  the  spherical  model  without  cracking 
(Timoshenk  and  Goodier,  1970)  is 


ffs  PqOT  sin(7sr) 

/(^)  =  “  e^’^(cos(7sr)  +  — ^ 


), 


where  r  =  t  -  (r  -  6)10^. 

The  displacement  history  at  r  (Timoshenk  and  Goodier,  1970)  is 
u{t,r)  =  3  +  TT— ((7S  +  - - )sin(7sr)  - 


C.  r 


s  rs 


), 


(3.13) 


(3.14) 


where  A  =  poa?/{2pyCc). 

The  first  term  in  the  solution  corresponds  to  the  static  solution,  the  second  is 
dynamic  term.  The  typical  displacement  history  from  Eq.(3.14)  is  shown  in  figure 
3.14  in  which  the  experimental  data  from  Shot  SWl  also  is  shown  (see  Chapter  4). 
Ftom  this  figure,  we  can  see  that  the  solution  from  only  pure  pulse  spherical  wave 
propagation  can  explain  the  P-wave  profile  from  the  experiments  without  cracking 
but  not  explain  the  observed  results  as  shown  in  figures  3.11  and  3.12. 

The  reasons  for  the  difference  may  include: 


1.  The  pressure  history  in  cavity  is  not  a  single-step-like  function.  Because  of  the 
large  mismatch  of  impedances  between  gas  in  cavity  and  rocks,  shock  waves 
generated  by  explosions  in  cavity  generally  reverberate  several  times  before  the 
pressure  in  cavity  reaches  its  final  value.  This  means  that  the  pressure  history  in 
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cavity  may  have  a  multi-step  history.  This  also  can  be  seen  from  the  experimental 
data  as  shown  in  figure  3.11  and  3.12; 

2.  The  cracking  processes  may  play  a  very  important  role  in  the  P-wave  induced 
displacement  history  because  radial  cracking  generated  by  the  pressure  in  cavity 
and  P-wave  propagation  changes  the  boundary  conditions  for  the  field.  When 
pressure  in  cavities  is  larger  than  some  critical  value,  radial  cracks  are  generated 
from  inner  surfaces.  The  radial  cracks  will  make  the  tangential  stress  be  either 
zero  if  the  explosion  products  do  not  penetrate  into  cracks  or  the  same  pressure 
as  in  the  cavity  if  the  explosion  products  do  penetrate  into  cracks.  Whenever  the 
explosion  products  penetrate  into  cracks  or  not,  the  tangential  stress  in  rocks 
has  a  large  jump  during  cracking  because  the  tangential  stress  is  tensional  before 
cracking  and  compressional  after  cracking.  These  rapid  changes  give  rise  to  shear 
waves.  Radial  cracking  also  generates  P-waves  because  cracking  is  equivalent  to 
increase  cavity  radius  (a  crack  opening  does  not  apparently  change  the  pressure 
in  cavity). 

Displacement  induced  by  the  P-waves  generated  from  explosions  and  cracking 

Basing  on  the  analysis  above,  we  construct  a  model  considering  the  reverberation  of 
shock  waves  in  cavity  and  the  P-waves  generated  from  radial  cracking  process. 

1.  Pressure  history 

In  this  model,  we  assume  that  the  pressure  history  is  divided  into  three  time 
periods,  i.e.  reverberation  dominated  period,  cracking  dominated  period  and 
steady  period.  The  pressure  history  acted  on  the  cavity  wall  can  be  given  by 

p{t)  =  pi{t){H{t)  —  H{t  —  to))  +  Pii{t){H{t  —  to)  —  H{t  —  too))  -^PiiiH{t  —  too), 

(3.15) 

where  p/(t),  p//(t)  and  p//j(t)  are  the  pressure  histories  that  correspond  to  the 
three  periods,  respectively.  H(t)  is  step  function,  to  and  too  are  the  start  and 
stop  time  of  crack  propagation.  The  pressure  history  for  each  period  in  cavities 
is  estimated  approximately  as 
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(I)  Reverberation  dominated  period: 

Shock  waves  induced  by  detonation  of  an  explosive  reverberate  in  a  cavity,  the 
pressure  acted  on  cavity  wall  is  multi-step  and  radial  cracks  have  not  yet  prop¬ 
agated  in  rock. 

During  this  period,  it  is  assumed  that  the  pressure  acted  on  the  wall  of  cavity 
between  two  adjacent  reverberations  is  constant,  we  have 

N 

Pi{t)  =  '£piH{t-ti),  (3.16) 

t=i 

where  Pi  is  pressure  increment  after  each  reflection,  tj  is  the  reflection  time  for 
ith  step  and  N  is  reverberation  times. 

(II)  Cracking  dominated  period: 

In  general,  there  are  two  kinds  of  cracks  induced  by  the  explosions  in  cavity.  The 
first  is  the  nucleation  of  new  cracks  and  growth  of  existing  cracks  (micro-cracks) 
related  to  P-wave  propagation,  the  second  type  of  cracks  is  the  macro-cracks 
that  start  from  the  cavity  wall  and  propagate  into  rock.  Micro-cracks  appears 
just  after  P-wave  if  stress  related  with  P-wave  is  higher  than  the  local  tensional 
strength.  Macro-crack  propagation  is  controlled  by  fracture  dynamics.  In  this 
part,  only  the  P-wave  generated  from  macro-cracking  is  considered  because  P- 
waves  generated  from  micro-cracks  can  not  be  separated  from  the  initial  P-waves 
and  also  have  small  amplitude. 

When  pressure  in  cavities  is  larger  than  some  critical  value,  radial  cracks  are 
generated  from  inner  surfaces.  The  tangential  stress  in  rocks  has  a  large  change 
during  cracking  because  the  tangential  stress  is  tensional  before  cracking  and 
compressional  after  cracking.  Radial  cracking  generates  P-waves  because  crack¬ 
ing  is  equivalent  to  increase  pressure  in  a  cavity. 

There  are  two  ways  to  approximately  estimate  the  effect  of  cracking  on  displace¬ 
ments. 

(a)  The  first  is  to  solve  the  differential  equation  with  moving  boundary  as 
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(3.17) 


fit) + 27(i)/'(t) + 

a{t)  p 


where 


7(t)  = 


l  —  2u  Cc 
1-1/  a^{t)  ’ 


(3.18) 


a^it)  =  a+  f  C“(r)dr, 

•/to 


(3.19) 


where  a^{t)  is  radius  of  the  crack  tip  to  the  center  of  the  cavity,  to  is  the  time 
at  which  cracks  start  to  propagate,  C‘  is  crack  propagation  velocity.  Numerical 
solution  may  be  applicable  to  Eq.(3.17)  but  the  moving  boundary  with  fracture 
is  also  very  difficult  to  deal  with. 

(b)  To  obtain  an  approximate  pressure  profile,  we  assume  that  the  pressure  in  a 
cavity  changes  basing  on  that  the  total  force  acted  on  the  spherical  surface  with 
radius  a^{t)  is  equal  to  the  force  acted  on  the  surface  with  radius  a. 


p{t)a^  =  po{a  +  f  C^dr)^, 
Jto 


(3.20) 


where  po  is  the  pressure  in  cavity  just  before  cracks  start  to  propagate  and  we 
assume  that  the  radial  stress  at  the  spherical  surface  with  radius  d^{t)  is  equal 
to  Po.  This  gives  the  equivalent  pressure  in  cavity  as 


p{t)  =  po(l  +  /  —drf,  t  >=  to, 

•/to  ^ 


(3.21) 


(III)  Steady  period: 

After  cracking  stops,  the  pressure  in  cavity  is  assumed  to  be  constant 


Pm  =  Piiitoo)- 


(3.22) 


2.  Solutions  for  the  three  periods 
(a)  Solution  for  period  I 


FVom  Eq.(3.14),  the  displacement  induced  by  ith  reverberation  at  r  is 


,  .  Ai  .  e  7  C'ex  .  ,  V  CeCos(7srj), 

^)  =  7  ”  77^  sin(75r<) - - - ^),  (3.23) 


where  Af  =  p,a^/(2p7Ce),  Tj  =  t-tj  —  {r  —  a)/ Cg  and  t,-  is  the  ith  reverberation 
time. 


the  resultant  solution  is  just  the  superposition  of  the  solutions  for  all  N  rever¬ 
berations.  The  displacement  field  induced  by  reverberations  at  r  is 


»  =  1 


(3.24) 


(b)  Solution  for  period  II 

In  order  to  obtain  an  analytic  expression,  the  crack  velocity  is  assumed  to  be 
the  two  simplest  cases  as: 

(1)  Case  a:  crack  velocity  remains  constant  during  crack  propagation. 

(2)  Case  b:  crack  velocity  decreases  linearly  with  time,  i.e.. 


C^it)  =  Q(1  -  (3.25) 

where  t^  is  the  total  cracking  duration,  Cg  is  the  initial  crack  propagation  velocity 
and  to  is  the  time  for  cracks  to  start  to  propagate. 

For  Case  (a),  the  pressure  history  in  period  II  is 

p(f)=Po(l  +  — (t-to))^  (3.26) 

a 

Substituting  Eq.(3.26)  into  Eq.(3.10),  the  solution  is  found  to  be 


where 


/(r)  =  Xo  +  Xir  +  X^Ki{t)  +  X^K^ir), 


Xo  =  Ao{^+  3' 

72  y  (1  4.  52) 


+ 


P2(3-  s^) 
7^1  +  52)2^’ 


(3.27) 


(3.28) 
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+  .^(1+52))’ 

(3.29) 

V  ^0P2 

"^2=  ^2  > 

(3.30) 

V  A  fPo  2pi  .  2p2(s^-3)^ 

-^5  ^0[^2  ^3(1  +  52)  ^4(1  +  ^2)2  h 

(3.31) 

y  A  f  Po  .  Pi(«^-1)  2p2(3s-5-^) 

®  725  735(1  +  52)  74(1+52)2 

(3.32) 

'^°“p(1  +  52)’ 

(3.33) 

Ki  =  6“'’'’’ cos  (757), 

(3.34) 

Ki  =  e"*’'’’ sin  (757), 

(3.35) 

Pi  —  ^Po - > 

a 

(3.36) 

P2=P0(^? 

a 

(3.37) 

and 

r-a 

r  =  t  V 

(3.38) 

Substituting  Eq.(3.27)  into  Eq.(3.7),  the  displacement  history  in 

period  II  is 

^  =  (?r:  +  tt)  +  (777  + -2-)^  +  + 

C.r  r^  C.r  r*'  r-' 


+(■— ^•  _^-—  +  ^)Ki  +  (^  -  (3.39) 


ar 


C.r 


For  Case  b,  the  pressure  history  is 

p  =  p<,(l  +  S((t-i„)-^i^))=.  (3.40) 

Substituting  Eq.(3.40)  into  Eq.(3.10),  the  solution  is 


/(r)  =  Co  +  ajT  +  a2T^  +  +  a4T^  +  a^Ki  +  a6-K^2,  (3-41) 


where 


33 


do  =  Ao{biPo  +  2^162  +  P2i>3(3  —  s^)  +  24p364(l  —  s^)  + 


+24^465(1  —  6s^  +  s^  +  4s  —  4s^)s  ^),  (3.42) 

oi  =  Ao(pibi  +  4p2b2  +  663(3  -  5^)  +  96^464(1  -  s^)),  (3.43) 

02  =  .^0(^261  +  6^362  +  12^463(3  —  5^)),  (3'44) 

fls  =  Aoipabi  +  8^462),  (3.45) 

04  =  ^oP46i,  (3.46) 


05  =  j4o(— 6ipo"l — 2^162  +  2(s^  —  3)p263  ~ 

24p3(l  —  5^)64  —  24p4(l  —  6s^  —  4s^  +  +  4s)bss~^),  (3.47) 


and 


06  =  Ao(-biPos  ^b2Pi(s^  -  l)s  ^  +  2b3P2(3s-s  ^)  -  664P3(s  ^-6s  +  s^)- 
-2465S-V4((1  +  s^)®  -  (1  +  4s  -  6s^  -  4s®  +  sy)°-^),  (3.48) 


where 


61  = 


(3.49) 


62  = 


1 

7^(1 +  s2)’ 


(3.50) 
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1 


®  7^(1  + S2)2’ 

(3.51) 

7®(1  +  52)3’ 

(3.52) 

65  -  ^ 

(3.53) 

«  2c§po 

Pi  =  , 

a 

(3.54) 

(3.55) 

(3.56) 

(3.57) 

Substituting  Eq.(3.41)  into  Eq.(3.7),  the  displacement  history  in  period  11  for 
Case  b  is 


2a2  Oi 


C,r  f^'  >2  ■  .  ^2. 


+ 


j  —  7<25  .  o-bs-r^  .  rOfi  'Vttfi  +  'vsaK. 

+  <  c,r  +  +  (;:?-  — cj  (3.58) 


(c)  Solution  for  period  III 

Substituting  Eq.(3.22)  into  Eq.(3.10),  the  solution  is 


~  ~  ®  ^‘(cos(7si)  - 


sin(7st)^ 


(3.59) 
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then,  substituting  Eq.(3.59)  into  Eq.(3.7),  the  displacement  is 


u^"{r) 

where 


Pia^  ,Ce  ,  -yr'rf'Y  ,  C*  COS(7St') 

— +  e  ^  ((-  +  7s- — )sm(7ST) - - - )).  (3.60) 

ZpyUeT  r  s  r$  r 


t'  =  T  -{■  till. 


(3.61) 


3.  Displacement  history  induced  by  P-waves 

Then  the  displacement  history  induced  by  P-waves  is 


u  =  -  H{t  -  fo))  +  u"{H{r  -  to)  -  H{t  -  too))  +  -  too),  (3.62) 


3.3.4  Comparison  between  the  observation  and  theoretical  predictions  of 
displacement  history  induced  by  P-waves 

Parameters  estimation 

Because  the  ratio  of  cavity  radius  to  explosive  radius  is  about  3,  there  is  no  separation 
between  shock  wave  in  air  in  the  cavity  and  the  adiabatic  release  wave  of  explosion 
products.  The  shock  wave  velocity  in  air  can  be  estimated  using  CJ  theory.  The 
particle  velocity  after  the  CJ  detonation  is 

Vp  =  (3.63) 

where  D  is  CJ  wave  velocity. 

For  PETN,  the  typical  CJ  wave  velocity  is  8  km/s  (Dobratz,  1974).  The  particle 
velocity  then  is  about  3.3  km/s.  The  vibration  period  in  the  cavity  is 

t„  R5  a/(2  ♦  3.3)  w  2. 3/^5.  (3.64) 

Prom  experimental  data,  the  rise-time  of  the  particle  displacement  induced  by  P- 
wave  is  ~  10  ps,  about  three  time  vibrations.  This  also  can  be  verified  from  the 
profiles  of  the  displacement  rise-front  from  the  experiments.  Then,  we  have  N=3, 
ti=0.0,  t2=4.6/xs,  t3=9.2  ^s,  Pi=P2=P3=Po/3.  Also,  to=10.0/xs. 

It  theoretically  is  found  that  the  maximum  velocity  of  a  mode-I  crack  is  Rayleigh 
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Figure  3.15:  Pressure  history 

wave  velocity  (Freund  and  Clifton,  1972).  Compared  to  the  crack  velocity  in  rocks 
from  seismology  (Hiroo  Kanomori,  1993),  the  fault  propagation  velocity  is  ~  70%  of 
the  theoretical  maximum  velocity  of  a  mode-II  crack,  i.e.,  shear  wave  velocity.  Crack 
propagation  velocity,  C*^,  is  assumed  to  be  70%  of  Rayleigh  wave  velocity,  C*,  because 
the  cracks  generated  in  the  experiments  are  mode-I.  Then,  too  is 


for  Case  a 

(3.65) 

2*  Lc 

“  0.7C« 

for  Case  b 

(3.66) 

where  Lc  is  crack  length. 

For  shot  SW4,  the  crack  length  with  the  explosion  product  trace  on  it  is  roughly 
about  50  mm,  so  too  38  //s  for  Case  a  and  65  fis  for  Case  b. 

Pressure  history  in  cavity 

Based  on  the  parameters  estimated  above,  the  pressure  profiles  normalized  by  po  for 
Shot  SW4  are  schematically  shown  in  figure  3.14. 

Displacement  history 

Figure  3.15  shows  the  calculated  results  of  displacement  time  history  based  on  the 
pressure  history  Case  a  and  Case  b.  From  this  figure,  it  seems  that  the  displacement 
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Figure  3.16:  Comparison  between  the  experimental  data  and  the  results  from  the  equivalent  models 
A:  Shot  SW4,  Gauge  2, 

B:  Constant  crack  velocity  model  (Eqs.(3.24)  and  (3.39)), 

C:  Constant  deceleration  model  (Eqs.(3.24)  and  (3.58)), 

D:  Constant  step  pressure  model  (Eq.(3.14)), 

E:  Show  SW9,  Gauge  4 


38 


history  recorded  from  experiments  is  between  the  two  extreme  conditions,  that  may 
reflect  the  real  situation  that  cracks  propagate  neither  at  constant  velocity  nor  at 
constant  deceleration.  Prom  the  comparison  between  the  theoretical  estimation  based 
on  Cases  a  and  b  and  the  experimental  data,  cracking  process  has  very  large  effect  on 
the  history  of  displacement  or  strain.  First,  cracking  processes  increase  the  amplitude 
of  low  frequency  component  because  cracking  is  equivalent  to  increase  the  pressure 
duration  in  a  cavity.  Second,  cracking  processes  increases  the  complexity  of  the 
directional  dependence  of  displacement  history  because  of  nonuniform  distribution  of 
cracks. 


3.3.5  P-wave  decay 

Prom  the  solutions  from  the  model  we  constructed  for  Period  I  as  given  in  Eq.(3.24), 
the  peak  displacement  at  r  induced  by  P-wave  is 


^max  — 


(3.67) 


where 

""  2^  ~  exp(-7Ti)(— +  cos(7STi)),  (3.68) 

and 

+  -)  sin(7STi)  exp(-7Ti).  (3.69) 

Ai  and  B\  can  be  determined  experimentally  using  the  measured  maximum  data 
at  different  distance  from  the  source.  The  best  fit  to  the  measured  data  is  shown  in 
figure  3.16.  The  best  fit  expressions  are 


Uraax  =  1.005  X  10V“2  -I-  676.6r"\  for  51^4,  (3.70) 


Umax  =  0.953  X  10V~2  530.1r~\  for  SW9,  (3.71) 

where  Umax  is  in  /im  and  r  is  in  mm. 

Based  on  the  relation  between  maximum  P-wave  particle  displacement  and  the 
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propagation  distance,  the  maximum  radial  strain,  e,.,  and  tangential  strain,  £<,  for 
SW4  are 

Br  =  =  -201r-"  -  0.677r-^  (3.72) 

=  100.5r-3  -f  0.677r-2. 
r 

(3.73) 

and  for  SW9 

£r  =  =  -190.6r-3  -  0.53r-2, 

(3.74) 

=  95.3r-3  +  O.SSr-^. 
r 

(3.75) 

Prom  Hooke’s  law,  we  have 

(3.76) 

(1  d-  i/)(l  -  2u) 

where  and  Ot  are  radial  and  tangential  stresses,  respectively.  Here,  v  is  Poisson’s 
ratio.  Based  on  the  experimental  data,  v  is  0.258  for  the  limestone  and  E  is  47.7 
GPa.  We  have  the  maximum  stresses, 

for  SW4: 

=  -1.19  X  10^(636r-®  -f  r-2), 

(3.78) 

at  =  3.92  X  10^(96.8r-^  +  r'^). 

(3.79) 

for  SW9: 

a,  =  -0.94  X  10‘‘(768r-2  d-  r'^), 

(3.80) 

at  =  3.06  X  10^(117r-^  d-  r'^). 

(3.81) 

3.4  Analysis  of  the  S-waves  Generated  from  Explosions  and 
Cracking 

3.4.1  S-wave  arrival 

Basically,  S-waves  are  generated  from  at  least  three  different  sources: 
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Figure  3.17:  Peak  P-wave  induced  particle  displacement  versus  propagation  distance 

1.  Micro-cracking. 

2.  shock  wave-to-S  conversion  at  cavity  wall. 

3.  Macro-cracking. 

In  the  case  of  each  of  these  three  sources,  the  arrival  time  of  the  S-wave  from  the 
second  source  can  be  estimated  with  relatively  high  accuracy  because  we  know  the 
P-wave  arrival  time  at  the  cavity  wall.  The  arrival  time  at  radial  distance,  r,  of  the 
S-wave  from  the  second  source,  Ts2,  is 


Ts2  = 


Ce-C, 

CgCa 


(r  -  a)  + 


(3.82) 


where,  C,  is  S-wave  velocity. 

The  duration  of  the  S-wave,  Ars2,  should  approximately  equate  to  the  rise  time 
of  P-waves. 

The  arrival  time  at  radial  distance,  r,  of  the  S-wave  from  the  first  source  is 

Tsl  =  -Rc)  +  (3.83) 
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where  Rc  is  the  radius  beyond  which  there  is  no  micro-cracks  induced  by  P-waves. 
The  duration  of  this  S-wave,  ATsl,  is 

^Tsl  =  -  <■)■  (3.84) 

The  accuracy  of  this  arrival  time  is  dependent  on  how  to  determine  /2c(in  the 
following  estimation,  is  chosen  to  be  50  mm  beised  on  the  experimental  data). 
The  arrival  time,  TsS,  of  the  S-wave  from  the  third  source  is 

^  (3.85) 

where  9c  is  the  angle  between  the  two  radial  lines  one  of  which  points  to  the  nearest 
macro-crack  generated  and  the  other  points  to  the  gauge  location.  The  last  term  on 
the  right  hand  side  in  Eq.  (3.85)  is  from  the  fact  that  the  S-wave  from  the  cracks 
that  is  not  along  the  line  connecting  the  gauge  with  the  center  of  the  cavity  arrives 
later.  The  earliest  and  last  arrival  times,  TsZ^  and  TsZi,  are  correspond  to  =  0 
and  6c  =  90®,  respectively.  Then,  we  have. 


r53,=  ^^(r-a)+(o, 

Ts3,  =  -  a)  + 1„  +  ~  (’■  -  “) 


(3.86) 


(3.87) 


CcC,  "  C. 

Using  the  expressions  above,  the  three  S-wave  arrivals  are  marked  on  the  records 
shown  in  figures  3.7  and  3.8.  Prom  these  figures,  the  predictions  are  in  very  good 
agreements  with  the  experiment  data  for  Ts2  and  Ts3.  The  data  clearly  demonstrates 
the  S-wave  generated  from  shock  wave-to-S  conversion  on  cavity  wall  and  cracking. 


3.4.2  S-wave  anaplitudes 

Figure  3.18  shows  the  maximum  strains  from  the  gauges  along  direction  2  (Fig.  2.1) 
for  shear  waves  generated  from  shock  wave-to-S  conversion  and  macro-cracking  at 
different  locations  for  Shot  SW4  and  SW9.  Due  to  the  difficulty  to  determine  the 
S-wave  incident  angle  at  each  gauges  directly,  the  two  gauges  at  the  same  location  is 
used  to  determine  the  incident  angle  approximately  and  then  the  peak  displacement 


42 


Distance  from  cavity  center  (mm) 


Figure  3.18:  Maximum  shear  wave  strain  versus  propagation  distance 

of  the  S-waves  are  determined.  The  procedure  to  determine  S-wave  incident  angle 
using  the  two  gauges  at  the  same  location  is:  (1)  to  find  out  the  ratio  of  the  strain 
from  the  gauge  along  direction  2  to  the  strain  from  the  gauge  along  direction  1;  (2) 
to  calculate  the  incident  angle  using  the  ratio  and  Eqs.(2.42)  and  (2.40).  We  must 
mention  that  this  method  is  not  very  precise  because  P-waves  interfere  with  S-waves 
but  at  least  it  gives  the  first-order  estimate  of  S-wave  amplitudes  induced  by  shock 
wave-to-S  conversion  at  cavity  wall  and  crack  propagation  in  rock.  Figure  3.19  shows 
the  maximum  displacement  of  S-waves  at  different  distance  from  source.  Prom  this 
figure,  it  seems  that  the  shear  wave  generated  from  shock  wave-to-S  conversion  has 
higher  regularity  than  the  shear  wave  generated  from  macro-cracking  regarding  the 
amplitude  decay  with  distance,  this  just  clearly  shows  the  complexity  of  shear  waves 
due  to  crack  propagation.  The  increase  of  the  S-wave  amplitude  generated  by  shock 
wave-to-S  conversion  at  cavity  wall  is  believed  to  be  from  the  source. 

3.4.3  S-wave  profiles  generated  from  cracking 

The  S-wave  amplitudes  and  profiles  recorded  by  different  gauges  are  dependent  on 
the  relative  orientation  of  the  gauges  with  respect  to  specific  cracks  (similar,  in  some 
respects,  to  the  directivity  of  wave  profile  in  seismology  ).  It  is  not  surprised  to 
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Figure  3.19:  Maximum  shear  wave  displacement  versus  propagation  distance 

see  that  the  S-wave  profiles  generated  from  cracking  vary  so  much  from  gauge  to 
gauge.  In  order  to  calculate  S-wave  profiles  generated  by  cracking,  we  assume  that  the 
shear  strain  decays  with  radial  distance  following  the  same  rule  as  P-waves  does(  this 
assumption  results  in  lower  estimation  of  S-wave  amplitude  because  P-wave  generally 
decays  much  faster  than  S-wave  does(Lay  and  Wallance,  1995)  but  does  not  have 
effects  on  the  wave  profiles,  and  we  only  consider  the  wave  profile  here),  then 

=  >  (3.88) 

i  ^ 

where  e,®  is  the  shear  strain  increase  generated  from  the  ith  crack  extending  to  radius 
a^ct(t),  Es  is  the  shear  strain  observed  at  r. 

In  order  to  obtain  an  order-of-magnitude  estimate  of  shear  waves  profiles  generated 
from  cracks,  we  assume  that  resultant  shear  strain  from  multi-cracks  is  equivalent  to 
that  from  a  single  crack  that  propagates  toward  to  gauges.  Then,  the  shear  strain 
profile  is  estimated  assuming  an  equivalent  crack  propagation  velocity.  We  still  only 
consider  the  two  extreme  cases,  i.e.,  crack  propagation  at  constant  velocity  and  at 
constant  deceleration. 
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For  constant  velocity  model,  the  shear  strain  profile  will  be 


e,  =  £!j((o  +  c;((  -  «.))/£)=■»,  (3.89) 

where  Cg  is  crack  propagation  velocity,  a  is  the  initial  cavity  radius  and  to  is  the 
initiation  time  of  cracks. 

For  constant  deceleration  model,  the  shear  strain  history  is 

£.  =  4((»  +  c;((  -  to)*" ~  *)/r)^\  (3.90) 

u 

where  t^  is  time  duration  of  crack  propagation. 

Based  on  the  two  equations  above,  the  shear  strain  history  is  calculated  as  shown 
in  figure  3.20  assuming  that  crack  propagation  time  is  same  for  both  model  in  order 
to  compare  them  together.  In  this  figure,  it  shows  that  experimental  results  from 
gauge  3  for  SW4  and  SW9  and  also  the  results  from  the  model  with  the  first  half  part 
at  constant  velocity  and  the  second  half  part  at  constant  deceleration.  It  seems  that 
a  constant  deceleration  or  a  constant  velocity  followed  by  a  constant  deceleration  fits 
the  data  better  than  the  constant  velocity. 

3.5  Conclusions 

Basing  the  experimental  results  and  the  analysis  above,  we  concluded 

1.  A  device  is  developed  to  drill  spherical  cavity  with  the  diameter  up  to  10  cen¬ 
timeters  in  rocks. 

2.  The  records  from  the  experiments  clearly  demonstrate  the  P-waves  generated 
from  both  explosions  and  crack  propagation.  The  peak  particle  displacement 
induced  by  P-waves  ranges  from  9  to  5  /rm  at  150  to  210  mm  from  source  for 
the  spherical  cavity  experiment  with  1.9  gram  explosives,  and  8  to  4  /im  at  150 
to  220  mm  from  source  for  spherical  experiment  with  0.8  gram  explosives.  The 
step-like  pressure  model  can  be  used  to  explain  the  experimental  results  on  the 
peak  P-wave  induced  particle  displacement  versus  propagation  distance  from  the 
explosions  in  spherical  cavities. 
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3.  More  important,  the  experimental  results  from  the  explosions  in  sphericcil  cavi¬ 
ties  verify  the  S-waves  generated  from  both  shock  wave-to-S  conversion  at  cavity 
wall  and  crack  propagation  in  rocks.  The  amplitude  of  S-waves  from  shock 
wave-to-S  conversion  is  from  10  to  27  /xm  at  150  to  180  mm  from  source  for  the 
spherical  experiment  with  1.9  gram  explosives,  8  to  21  ^m  at  150  to  180  mm 
from  source  for  the  spherical  experiment  with  0.8  gram  explosives.  The  ampli¬ 
tude  of  S-waves  generated  by  crack  propagation  is  from  9  to  19  ^m  at  150  to 
180  mm  from  source  for  the  spherical  experiment  with  1.9  gram  explosives,  8  to 
19  /xm  at  150  to  180  mm  from  source  for  the  spherical  experiment  with  0.8  gram 
explosives.  The  ratio  of  S-wave  to  P-wave  amplitudes  is  about  2. 

4.  A  theoretical  model  including  multi-reverberation  and  crack  propagation  is  pro¬ 
posed  to  explain  the  differences  on  P-wave  induced  displacement  history  between 
the  observed  and  the  step-like  pressure  model.  The  predicated  profiles  from  this 
model  fits  the  experimental  data  qualitatively  as  shown  in  Fig.  3.16. 

5.  A  simple  dynamical  tensional  strength  model  is  used  to  estimate  the  S-wave 
profile  generated  from  crack  propagation  in  rock.  The  S-wave  profile  predicted 
from  this  model  fits  the  experimental  data  quantitatively. 
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Chapter  4  Wave  Generation  From  Explosions  in 
Cylindrical  Cavities  in  Rock 


4.1  Introduction 

Initially  to  test  the  apparatus  constructed  and  developed  in  Chapter  1,  we  conducted 
two  experiments  using  a  simple  cylindrical  charge  detonated  in  a  cylindrical  hole 
centered  in  the  rock  sample  depicted  in  figures  3.3  and  4.1.  At  same  time,  the  cavity 
geometrical  effects  on  P-  and  S-waves  generated  from  explosions  can  be  investigated 
from  the  experimental  data  on  spherical  and  cylindrical  cavities. 

4.2  Experiments  and  Results 

4.2.1  Experimental  Detadls 

Two  cylindrical  test  experiments  were  conducted  on  Bedford  limestone.  One  was  a 
decoupled  explosion,  SWl,  and  the  other  a  tamped  explosion,  SW2.  The  parameters 
used  in  the  experiments  and  the  equilibrium  pressure  calculated  using  Eq.(3.1)  are 
listed  in  Table  4.1. 

The  rock  sample  is  assembled  using  two  blocks  as  shown  in  figure  4.1.  The  rock 
sample  with  strain  gauges  is  placed  inside  a  tank  pressurized  to  10  bar  as  shown  in 
figure  3.3.  To  prevent  the  sample  from  becoming  invaded  by  water,  a  plastic  film  [ 
0.21  mm  thickness  (Vinyl  film  sheet  No.  8562K5,  Warp  Bros.)]  is  wrapped  on  the 
surface  of  the  rocks.  The  gauges  and  explosives  used  in  experiments  are  identical  to 
that  used  in  spherical  experiments. 

4.2.2  Experimental  Results 

The  recorded  strains  for  the  two  experiments  are  shown  in  figures  4.2  for  Shot  SWl 
and  figure  4.3  for  Shot  SW2.  From  the  recovered  rock  sample,  there  is  minimal  damage 
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Strain  gauges  (7,  8) 


5.6 

3,4 

1,2 

Rock  joint 

Explosives 

Rock  sample 


(A)  Layout  of  rock  sample,  strain  gauges  and  explosive  charge 


(B)  Sketch  of  polarization  directions  of  strain  gauges 


Figure  4.1:  Sample  and  strain  gauge  layout  for  cylindrical  cavity  experiments 


49 


Table  4.1:  Data  for  cylindrical  cavity  experiments 


Shot  # 

Dimensions 
of  explosives 

Mass  of 
explosives 

Type  of  explosives 

Dimensions 
of  cavities 

SWl 

4>  5x6  mm 

0.24  g 

Pellet  pressed  PETN 

030x30  mm 

55 

(j)  10x10  mm 

1.2  g 

Deta  sheet  PETN 

0  12x12  mm 

4420 

Figure  4.2:  Experimental  results  of  the  decoupled  explosion,  Shot  SWl 
The  signals  of  gauge  7  and  8  are  very  noisy  and  are  not  shown  here. 

near  the  cavity  wall  in  the  target  used  for  the  decoupled  experiment,  and  there  are 
many  of  the  cracks  with  varying  length  up  to  122  mm  and  different  orientations  in 
the  target  used  for  the  tamped  experiment. 

Prom  figure  4.2,  it  can  be  seen  that  the  wave  profiles  from  the  decoupled  experi¬ 
ment,  SWl,  is  very  similar  to  the  results  from  the  analytic  solution  shown  in  figure 
3.13  although  the  cavity  is  not  spherical.  Prom  figure  4.3,  the  wave  profiles  from 
the  tamped  experiment,  SW2,  have  the  typical  characteristics  of  the  wave  profiles 
from  the  spherical  experiments.  The  difference  between  the  decoupled  and  tamped 
experiments  once  again  demonstrates  the  role  of  cracking  in  the  generation  of  P-  and 
S-waves  because  there  is  minimal  damage  in  the  decoupled  experiment  and  the  rock 
sample  is  heavily  cracked  in  the  tamped  experiment. 


Figure  4.3:  Experimental  results  of  the  tamped  explosion,  Shot  SW2 


Figure  4.4:  P-wave  amplitudes  from  decoupled  experiment 

4.3  P-waves  Generated  from  Explosions 

4.3.1  P-waves  generated  from  the  decoupled  experiment 

Figure  4.4  gives  experimental  results  on  the  P-wave  profiles  generated  from  the  decou¬ 
pled  experiment.  The  peak  particle  displacements  from  P-wave  are  shown  in  figure 
4.6.  As  mentioned  above,  although  the  P-wave  profile  shape  can  be  explained  quan¬ 
titatively  using  the  analytic  solution  based  on  spherical  source  assumptions  as  shown 
in  figure  3.13,  experimental  results  show  the  direction-dependent  of  the  wave  shape. 
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Figure  4.6:  P-wave  amplitude  versus  propagation  distance  for  SWl  and  SW2 

4.4.2  S-waves  generated  from  the  tamped  experiment 

Figure  4.8  gives  the  arrivals  of  the  three  S-waves  from  the  three  different  source  as 
discussed  in  section  3.4.  As  the  same  as  in  the  spherical  experiments,  we  assume 
that  the  macro-cracks  start  to  propagate  when  P-wave  reaches  the  peak  value.  The 
micro-crack  radius  is  chosen  to  be  122  mm  based  on  the  cracks  in  the  recovered 
target.  Although  the  source  is  cylindrical,  there  still  is  a  good  agreement  between 
the  predicted  arrivals  of  the  S-waves  and  the  experimental  results.  We  can  see  the 
S-waves  from  the  shock  wave-to-S  conversion  and  macrocracks  on  all  the  records. 
Figure  3.9  gives  the  maximum  strain  related  with  the  shear  waves  generated  from 
shock  wave-to-S  conversion  and  macrocracking  along  direction  2. 


4.5  Conclusions 

ifrom  the  experimental  results  of  the  tamped  and  decoupled  explosions  in  cylindrical 
cavities,  we  concluded 

1.  The  records  verified  that  the  characteristics  of  the  signals  from  the  gauges  at 
different  locations  predicted  in  Chapter  1.  It  means  that  method  developed  in 
Chapter  1  can  be  used  to  monitor  P  and  S-wave  generations  from  explosions  in 
cavities  in  rock. 
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Figure  4.6:  P-wave  amplitude  versus  propagation  distance  for  SWl  and  SW2 
4.4.2  S-waves  generated  from  the  tamped  experiment 

Figure  4.8  gives  the  arrivals  of  the  three  S-waves  from  the  three  different  source  as 
discussed  in  section  3.4.  As  the  same  as  in  the  spherical  experiments,  we  assume 
that  the  macro-cracks  start  to  propagate  when  P-wave  reaches  the  peak  value.  The 
micro-crack  radius  is  chosen  to  be  122  mm  based  on  the  cracks  in  the  recovered 
target.  Although  the  source  is  cylindrical,  there  still  is  a  good  agreement  between 
the  predicted  arrivals  of  the  S-waves  and  the  experimental  results.  We  can  see  the 
S-waves  from  the  shock  wave-to-S  conversion  and  macrocracks  on  all  the  records. 
Figure  3.9  gives  the  maximum  strain  related  with  the  shear  waves  generated  from 
shock  wave-to-S  conversion  and  macrocracking  along  direction  2, 


4.5  Conclusions 

From  the  experimental  results  of  the  tamped  and  decoupled  explosions  in  cylindrical 
cavities,  we  concluded 

1.  The  records  verified  that  the  characteristics  of  the  signals  from  the  gauges  at 
different  locations  predicted  in  Chapter  1.  It  means  that  method  developed  in 
Chapter  1  can  be  used  to  monitor  P  and  S-wave  generations  from  explosions  in 
cavities  in  rock. 
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Figure  4.9:  S-wave  amplitude  versus  propagation  distance  for  SWl  and  SW2 

2.  The  peak  particle  displacement  induced  by  P- waves  ranges  from  35  to  10  ^m  at 
150  to  240  mm  from  source  for  the  tamped  cylindrical  cavity  experiment  with  1.2 
gram  explosives  and  7  to  2  /tm  at  150  to  190  mm  from  source  for  the  decoupled 
cylindrical  cavity  experiment  with  0.24  gram  explosives.  The  step-like  pressure 
model  can  not  be  used  to  explain  the  experimental  results  on  the  peak  P-wave 
induced  particle  displacement  versus  propagation  distance  from  the  explosions 
in  cylindrical  cavities.  This  results  demonstrates  the  cavity  geometrical  effects 
on  wave  pattern. 

3.  The  shear  waves  generated  from  shock  wave-to-S  conversions  at  cavity  wall  were 
observed  to  be  between  50  and  400  along  direction  2  at  150  to  250  mm  from 
source  for  tamped  explosions.  The  shear  waves  induced  by  crack  propagation 
were  recorded  to  be  between  100  and  250  fiE  along  direction  2  at  150  to  250  mm 
from  source  for  the  tamped  experiment. 
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Chapter  5  Summary 


Basing  on  the  methods  developed  and  the  experimental  results  and  analysis  given 
above,  we  briefly  conclude: 

1.  A  new  measuremental  method  to  monitor  the  P-  and  S- waves  with  very  high 
frequency  has  been  developed.  The  experimental  results  verify  the  predicted 
strain  characteristics  of  P-  and  S-  waves  generated  from  explosions  and  cracking. 

2.  A  device  is  developed  to  drill  spherical  cavity  with  the  diameter  up  to  10  cen¬ 
timeters  in  rocks  like  limestone. 

3.  The  records  from  the  experiments  clearly  demonstrate  the  P- waves  generated 
from  both  explosions  and  crack  propagation.  The  peak  particle  displacement 
induced  by  P-waves  ranges  from  9  to  5  fim  at  150  to  210  mm  from  source  for 
the  spherical  cavity  experiment  with  1.9  gram  explosives,  and  8  to  4  /zm  at 
150  to  220  mm  from  source  for  spherical  experiment  with  0.8  gram  explosives, 
35  to  10  /xm  at  150  to  240  mm  from  source  for  the  tamped  cylindrical  cavity 
experiment  with  1.2  gram  explosives  and  7  to  2  /zm  at  150  to  190  mm  from  source 
for  the  decoupled  cylindrical  cavity  experiment  with  0.24  gram  explosives.  The 
step-like  pressure  model  can  be  used  to  explain  the  experimental  results  on  the 
peak  P-wave  induced  particle  displacement  versus  propagation  distance  from 
the  explosions  in  spherical  cavities  but  not  in  cylindrical  cavities.  This  results 
demonstrates  the  cavity  geometrical  effects  on  wave  pattern.  The  experimental 
results  from  both  spherical  and  cylindrical  cavities  show  that  P-wave  profiles  are 
similcir  with  each  other. 

4.  More  important,  the  experimental  results  in  both  spherical  and  cylindrical  cavi¬ 
ties  verify  the  S-waves  generated  from  both  shock  wave-to-S  conversion  at  cavity 
wall  and  crack  propagation  in  rocks.  The  amplitude  of  S-waves  from  shock 
wave-to-S  conversion  is  from  10  to  27  /zm  at  150  to  180  mm  from  source  for  the 
spherical  experiment  with  1.9  gram  explosives,  8  to  21  /zm  at  150  to  180  mm 
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from  source  for  the  spherical  experiment  with  0.8  gram  explosives.  The  ampli¬ 
tude  of  S-waves  generated  by  crack  propagation  is  from  9  to  19  fim  at  150  to 
180  mm  from  source  for  the  spherical  experiment  with  1.9  gram  explosives,  8  to 
19  /xm  at  150  to  180  mm  from  source  for  the  spherical  experiment  with  0.8  gram 
explosives. 

5.  A  theoretical  model  (Section  3.3.4)  including  multi-reverberation  and  crack  prop¬ 
agation  is  proposed  to  explain  the  differences  on  P-wave  induced  displacement 
history  between  the  observed  and  the  step-like  pressure  model.  The  P-wave  pro¬ 
files  predicted  from  this  model  fits  the  observation  qualitatively  as  shown  in  Fig. 
3.16. 

Prom  the  results  presented  above,  we  believe  that  the  further  study  of  the  charac¬ 
teristics  of  P-  and  S-wave  profiles  generated  from  explosions,  shock  wave-to-S  conver¬ 
sion  and  crack  propagation  driven  by  explosive  products  may  lead  to  the  development 
of  a  criterion  to  discriminate  underground  chemical  explosions  from  underground  nu¬ 
clear  explosions  because 

1.  cracks  generated  in  chemical  explosions  may  be  much  longer  than  that  in  nuclear 
explosions  as  shown  in  Fig.  5.1.  Longer  crack  length  will  result  in  both  P-  and 
S-waves  with  more  lower  frequency  components. 

2.  the  conversion  of  shock  wave-to-S  wave  at  cavity  wall  may  only  happen  in  chem¬ 
ical  explosions. 

3.  P-wave  profiles  contain  the  informations  about  cavity  size. 
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Chemical  explosion  Nuclear  explosion 


Figure  5.1:  Possible  damage  difference  between  chemical  and  nuclear  explosions 
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